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There exists a variety of physically interesting situations described by continuous maps that are 
nondifferentiable on some surface in phase space. Such systems exhibit novel types of bifurcations 
in which multiple coexisting attractors can be created simultaneously. The striking feature of these 
bifurcations is that they lead to fundamentally unpredictable behavior of orbits as a system param- 
eter is varied slowly through its bifurcation value. This unpredictability gradually disappears as the 
speed of variation of the system parameter through the bifurcation is reduced to zero. 
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In the literature dealing with bifurcation theory, dif- 
ferentiable maps are commonly considered. On the other 
hand, maps that are piecewise differentiable occur in a 
variety of physical situations. We will consider here two 
dimensional maps. By a piecewise differentiable two di- 
mensional map we mean that the map is continuous and 
there is a curve (which we call the border) separating the 
phase space into two regions, such that the map is dif- 
ferentiable on both sides of the border, but not on it. 
In particular, in this letter we shall be concerned with 
the case where the map's derivative changes discontinu- 
ously across the border. This circumstance leads to a rich 
class of bifurcation phenomena called border- collision bi- 
furcations examples of which may be found in in- 
stances of grazing impact in mechanical oscillators ||-(7) , 
in piecewise-linear electronic circuits || , and in commer- 
cial power electronic devices |P,|Io|. 

In this letter, we study the occurrence and prevalence 
of border-collision bifurcations in which multiple coex- 
isting attractors are simultaneously created at the bi- 
furcation. We call these bifurcations multiple attractor 
bifurcations. 

Our principal point in this letter is that multiple at- 
tractor bifurcations display a striking new type of ex- 
treme sensitivity to noise. In particular, assume that 
the noise SN is at all times limited by some small value 
e > \SN\. Assume that before the bifurcation, the orbit 
follows a particular attractor (to within the noise) and 
then imagine that, as the system orbit evolves, the bifur- 
cation parameter is very slowly varied, eventually passing 
through its value at the bifurcation. The question is: Af- 
ter the bifurcation, which attractor does the orbit follow? 
We argue that no matter how small the noise level e may 
be, this question may be fundamentally unanswerable. 

As a simple representative example, we consider the 
electronic circuit with an analog switch pj , S, controlled 
by a comparator |12], O, as depicted in Fig. la. Let q 
denote the electric charge on the capacitor C, and let i 
denote the current flowing through the inductor L and 



the resistor R. For small values of the current, i, the 
voltage drop, iR, across the resistor R is smaller than 
the reference voltage, v re f. Consequently, the compara- 
tor, O, keeps the switch, S, open (off). However, when 
the current increases so that i > i cr n — v re f / R, the com- 
parator output changes state, thereby closing the switch, 
S. Ideally, the transition of the switch from an open 
state to a closed state (and vice versa), should be sharp 
and instantaneous. In reality, however, the situation is 
somewhat compromised by a variety of factors including 
finite gain and slew rate in the comparator op-amp, and 
finite rise and fall times of the switch. This causes a 
slight blurring of the i — i C rit boundary. Fortunately, the 
precision and speed afforded by present-day components 
restricts this boundary broadening to levels insignificant 
compared to the typical value of i in the circuit, thereby 
justifying the sharp boundary model. The entire sys- 
tem is driven by an externally applied voltage signal v(t) 
which is a square wave of amplitude v amp with a con- 
stant bias Vbias- The circuit model can be described by 
the following normalized equations: 



dQ/dT 
dl/dT 



I (I < 1) 

I- (Q + Qbias)/p (I>1) 

-n 2 Q-ri + F (o<r<i/2) 

-n 2 Q-TI-F (1/2<T<1) 



(1) 
(2) 



where / = i/i C rit, Q = (fq/ierit) - Qbias, and the time 
T is normalized to the period of the square wave drive, 
T = ft where / is the drive frequency. The parame- 
ters in (1) and (2) are p = fCR c , F = v amp (f 'Li crit Y x , 
f2 _1 = fVLC, r = R/(fL), and Q bias = v bias fC/i crlt . 
Figures 2(a) and (b) show bifurcation diagrams for this 
system for the two-dimensional time-1 map of this system 
(for parameters see the figure caption). The interesting 
fact is that when (as was done to get Figs. 2) one creates 
bifurcation diagrams by following the attractor as the 
bifurcation parameter is varied, upon addition of even 
extremely small noise (the level of which can be judged 
from the small thickening of the lines in the bifurcation 
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diagram), we find that after the bifurcation the orbit can 
go to different attractors: either a period three attractor 
(Fig. 2(a)) or a chaotic attractor (Fig. 2(b)). 

Before proceeding further, we introduce a canonical 
form Jl3| for border-collision bifurcations given in 
that allows us to treat piecewise-differentiable systems in 
a general system-independent manner. Let G M be a one- 
parameter family of piecewise smooth maps from the two 
dimensional phase space to itself, depending smoothly on 
the parameter /i. Let denote a fixed point (period one 
orbit) of defined on (j,q — a < fi < /io + a and depend- 
ing continuously on /i, for some a > (The analysis may 




where t{x) and d{x) are piecewise constant with a discon- 
tinuity at x = 0, that is, (t(x), d(x)) = (r<, d<) for x < 
and (t(x), d(x)) = (r > ,d>) for x > where r<,c? < ,r > , 
and d > are constants. In Eq. (3) the border has been 
transformed to the y-axis, /io = and = (0,0). We 
refer to (3) as the canonical form of the original map . 
By the piecewise linearity of (3) in (x, y) and its linearity 
in /i, it is invariant under (x,y,/j,) — > (Xx,Xy,Xfi). Thus, 
as \x is reduced towards zero, any attractor of the system 
and its basin of attraction must contract linearly with /i, 
with the attractor collapsing to the point (x,y) — (0,0) 
(cf., Figs. 2(c) and (d) in which both the period three at- 
tractor and the chaotic attractor collapse onto the point 
(0,0)). Thus it appears that, no matter how small the 
noise level e may be, there always exists a finite, positive 
/^-interval where the attractors are close enough for the 
noise to induce inter-attractor hopping, thereby render- 
ing prediction of the final attractor impossible @,^|| ■ 

The canonical form (3) depends on four parameters 
r < ,d < ,T > , and d> (and the bifurcation parameter ji). 
To find them from the original map G M , we evaluate its 
Jacobian matrix of partial derivatives for /i = /io on both 
sides of the border, infinitesimally close to E^ . The pa- 
rameters t < , t> and d < , c?> are then the trace and the 
determinant of the Jacobian matrix on the two sides of 
the border. As an example, we evaluate the canonical 
form for the system, Eqs. (1, 2), with parameters corre- 
sponding to those of Figs. 2(a) and 2(b). Figures 2(c) 
and 2(d) show bifurcation diagrams obtained from this 
canonical form (map G^) along with some very small ad- 
ditive noise. Note that this figure reproduces the partic- 
ular bifurcation phenomenon: there is a border-collision 
bifurcation of the fixed point attractor (/i < 0) to two 
attractors (/i > 0), one a period three attractor and the 
other a chaotic attractor. This is an example of a multi- 
ple attractor bifurcation. 

In our example, Eqs. (1, 2), the Jacobian determinant 
is constant, <i< = d > = d > 0, and for simplicity we 



be easily extended to include points on a period-p orbit 
by choosing the map to be the p-th iterate of G M ). As the 
parameter fi is increased from /io — a, we suppose that 
the fixed point E^ collides with the border at fx = /.to- 
We could also allow the border to depend on [i, but co- 
ordinates can always be chosen so that it remains fixed. 
For /i near /io the attractors are located near £^ , and 
one can expand the map to first order for (x, y) near E^ 
and (/z — /io) small. Nusse and Yorke [Q show that after 
suitable changes of coordinates and rescaling of the pa- 
rameter, the resulting continuous, piecewise linear map 
can be cast in the form 



henceforth limit our considerations to this case. Thus, 
for fixed d, the type of border-collision bifurcations that 
occurs depends on the two parameters r< and r> . We 
wish to explore this parameter space to assess the preva- 
lence of multiple attractor bifurcations. Figure 3 shows 
the (t<,t> )-plane for 0.5 < t< < 0.9, -2.2 < t> < -1.8, 
d = 0.3 and fi = 1. Two specific initial conditions, 
(x ,y ) = (7.0,2.0) and (x ,y ) = (1.2,0.0), were cho- 
sen for our simulations at every (t<,t>) grid point. The 
results show four distinct regions — two where both the 
initial conditions lead to a period-3 attractor, another 
where both the initial conditions lead to a chaotic attrac- 
tor, and still another where one of the initial conditions 
leads to a period-3 attractor while the other leads to a 
chaotic attractor. Clearly, the fourth region represents 
parameter values where at least two attractors co-exist 
after bifurcation, thus making multiple attractor bifur- 
cations possible. 

It has been shown so far how a novel form of unpre- 
dictability might arise in piecewise smooth maps. How- 
ever, this new uncertainty gradually gives way to cer- 
tainty as the bifurcation speed (the rate at which the 
bifurcation parameter is varied in time) is reduced to ar- 
bitrarily low values. In order to understand this consider 
the case where there are two attractors A and B, and fi 
and the upper bound e on the noise level are held fixed. 
For very small /z, the (noiseless) attractors are closer to- 
gether than the noise level e. For sufficiently large /i the 
attractors are so far apart that the noise cannot kick an 
orbit on one of the attractors out of its basin, and the 
orbit stays near the same attractor, A or B, forever. In 
general, however, there will be intermediate fj, values fj,A 
and hb such that hopping from A to B becomes impos- 
sible when /i > /i^, and hopping from B to A becomes 
impossible when ji > [Ag. Generically (in the absence of 
special symmetries) fiA 7^ M-Bj an( i we assume /is > /M; 
e.g., attractor B might be closer to the basin boundary 
than is attractor A. Therefore, if the bifurcation param- 
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eter were to be varied infinitely slowly, during the time 
when fiB > M > fJ-A, the system would certainly end up 
in A and it would remain there. Thus, in the quasistatic 
limit, there is no unpredictability as to which attractor 
the system ends up in. 

Thus unpredictability occurs at nonzero bifurcation 
speed and we now ask whether we can predict anything 
about the relative probabilities of ending up on the dif- 
ferent attractors given the bifurcation speed. 

For simplicity, let us consider the situation with just 
two post-bifurcation attractors, A and B, as described in 
a previous paragraph. Let v denote the bifurcation speed, 
while a(jx) and b([i) represent the probabilities of the sys- 
tem being in attractors A and B respectively. Then for 
small v we have 

a + b=l, (4) 
db/dt = -XsAb + Aab&, (5) 

where Xab(^) and Xba{ij) are the noise induced transi- 
tion probabilities from A to B and B to A, respectively. 
Since Xab{h) = for jj, > ha, we immediately have 

b[p) = b{^ A )e , (6) 

for fi > /ig, where we assume that we begin at some /i < 
li A- It should be pointed out that as v — > 0, - * 
as well, but only as v@ for Xab{h — * Ha) ~ ip-A — h)^ ■ 
Therefore, to leading order, > /is), the probability 
of the system ending up on the non-quasistatic attractor, 
satisfies ln[6(/i > ^b)] ~ —1/v. 

Two sets of simulations were conducted for Eq. (3) 
at the parameter values corresponding to the electronic 
circuit described before. The system was subjected to 
additive noise with bounded amplitude. In one case, the 
noise vector distribution had the form of a uniformly 
filled square centered at the origin. In the other case, 
the distribution was a uniformly filled circular disk also 
centered at the origin. In both cases, when the speed 
of variation of the bifurcation parameter was sufficiently 
small, the system tended to relax to the period three at- 
tractor. Figure 4 confirms that, at intermediate speeds, 
the probability of not going to the period three attractor 
indeed conforms to Eq. (6). 

In conclusion, in this paper, we have discussed a novel 
form of uncertainty that can arise in piecewise smooth 
systems. We have also pointed out how that uncertainty 
slowly fades into certainty as the speed of bifurcation is 
reduced to zero. 
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FIGURE CAPTIONS 
Figure 1 

(a) The Border Collision Circuit. 

(b) i < i crit . 

(c) i > icrit. 

Figure 2 

(a) and (b) Bifurcation diagrams for the two-dimensional time-1 map of the circuit depicted in Fig. 1 with Qbias = 

— 1.0, p = 0.10742, ft = 1.0642 and T = 1.2040. In creating the bifurcation diagrams (a) and (b) small noise was 
inserted leading to two different realizations. 

(c) and (d) Bifurcation diagrams for the equivalent canonical form, Eq. (3), near criticality with t< = 0.7, t> = —2.0 
and d = 0.3. In creating (c) and (d), small noise was inserted leading to two different realizations. 

Figure 3 

Two-dimensional period plot in parameter space for the canonical form, Eq. (3), with 0.5 < r< < 0.9, —2.2 < r> < 

— 1.8, d = 0.3, fi = 1, and initial conditions (x ,y ) = (7.0, 2.0) and (x ,y a ) = (1.2, 0.0). 

Figure 4 

Probability > ^2), plotted logarithmically, as a function of 1/v, based on 1,000,000 experiments for every value 
of 1/v. Solid lines and plus symbols represent square noise, while dashed lines and asterisks represent circular noise. 
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